Resonant backward scattering of light by a two-side-open subwavelength metallic slit 
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The backward scattering of TM-polarized light by a two-side-open subwavelength slit in a metal 
film is analyzed. We show that the reflection coefficient versus wavelength possesses a Fabry-Perot- 
like dependence that is similar to the anomalous behavior of transmission reported in the study 
[Y. Takakura, Phys. Rev. Lett. 86, 5601 (2001)]. The open slit totally reflects the light at the 
near-to-resonance wavelengths. In addition, we show that the interference of incident and resonantly 
backward-scattered light produces in the near-field diffraction zone a spatially localized wave whose 
intensity is 10-10 3 times greater than the incident wave, but one order of magnitude smaller than 
the intra-cavity intensity. The amplitude and phase of the resonant wave at the slit entrance and 
exit are different from that of a Fabry-Perot cavity. 
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agating modes inside the slit waveguide iSiifliiiiiSii^ At the 

resonant wavelengths, the transmitted field increases via 
the strong coupling of an incident wave with the waveg- 
uide modes giving a Fabry-Perot-like behavior ^iiiiiSi 
In the case of films, whose thickness are too small to sup- 
port the intra-cavity resonance, the extraordinary trans- 



I. INTRODUCTION 

The most impressive features of light scattering by 
subwavelength metallic nanostructures are resonant 
enhancement and localization of the light by excita- 
tion of electron waves in the metal (for example, see . . . . , , . , 
rcf: 1.2.3.4.5.6.7.8.9.10.11.12.13.14.15.16.17.18.19.20.21^ mechanism, the genera- 

T111 ,„ n r , t ^ tion of resonant surface plasmon polaritons and coupling 

In the last few years, a great number of studies have 

been devoted to the nanostructures in metal films, 

namely a single aperture, a grating of apertures and 

an aperture surrounded by grooves. Since the recent 

paper of Ebbesen and colleagues^ on the resonantly 

enhanced transmission of light observed for a 2D 

array of subwavelength holes in metal films, the 

resonant phenomenon is intensively discussed in the , . ,, ,, . „„ , . . .. , - .... 

21,22, 23.24.25,26,27,fem.gffMiag3tetiafe(j;§pectmg TE-polarized lights The 

i^imta^nasmonsEmcniayleigh anomalies were involved 



of them into radiationA^iiiiiSi 2 ^ Both physical mecha- 
nisms play important roles in the extraordinary transmis- 
sion through arrays of two-side-open slits (transmission 
gratings) and the resonant reflection by arrays of one- 
side-open slits (reflection gratings). A model of trapped 
(waveguide) modes has been recently used to show that 
an array of two-side-open slits can operate like a reflec- 



liter ature 4 - 5 - 6 - 7 - 8 - 9 - 1 0- 11 - 12 - 13 - 14 - 15 - 16 - 17 - 18 - 19 - 20 
Such a kind of light scattering is usually called a Wood's 
anomaly. In the early researches, Hessel and Oliner 
showed that the resonances come from coupling between 
nonhomogencous diffraction orders and eigenmodes of 
the grating^ Neviere and co-workers discovered two 
other possible origins of the resonances Al One appears 
when the surface plasmons of a metallic grating are 
excited. The other occurs when a metallic grating is 
covered by a dielectric layer, and corresponds to guided 
modes resonances in the dielectric film. The role of 
resonant Wood's anomalies and Fano's profiles in the 
resonant transmission were explained in the study. 8 

The phenomena involved in propagation through hole 
arrays are different from those connected with slit arrays. 
In a slit waveguide there is always a propagating mode in- 
side the channel, while in a hole waveguide all modes are 
evanescent for hole diameters smaller than approximately 
a wavelength. In the case of slit apertures in a thick metal 
film, the transmission exhibits enhancement due to a 
pure geometrical reason, the resonant excitation of prop- 



in explanation of reflective properties of such a kind of 
gratings 

The studiesiiiiiSiiS have pointed out that the ori- 
gin of anomalous scattering of light by a grating of 
slits (holes) can be better understood by clarifying the 
transmission and reflection properties of a single sub- 
wavelength slit. Along this direction, it was already 
demonstrated that the intensity of TM-polarized light 
resonantly transmitted through a single slit can be 10- 
10 3 times higher than the incident wavei*2Ai^ and that 
the transmission coefficient versus wavelength possesses 
a Fabry-Perot-like behavioniiii 4 ^. Unfortunately, the 
reflection properties of the slit have received no atten- 
tion in the literature. The very recent studjsii only con- 
cerned the problem by regarding the total reflection of 
TE-polarized light by a grating of two-side-open slits to 
properties of the independent slit emitters. 

In this article, the backward scattering of light by a 
two-side-open subwavelength slit is analyzed. To com- 
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pare properties of the light reflection with the extraor- 
dinary transmissioni^iai^, we consider the scattering of 
TM-polarized light by a slit in a thick metallic film of per- 
fect conductivity. From the latter metal property it fol- 
lows that surface plasmons do not exist in the film. Such 
a metal can be described by the Drude model for which 
the plasmon frequency tends towards infinity. The tradi- 
tional approach based on the Neerhoff and Mur solution 
of Maxwell's equations is used in the computations^^ 
The article is organized as follows. The theoretical back- 
ground, numerical analysis and discussion are presented 
in Section II. The summary and conclusions are given 
in Section III. The brief description of the model is pre- 
sented in the Appendix. 



II. NUMERICAL ANALYSIS AND DISCUSSION 

It is well known that when a light wave is scattered by 
a subwavelength metallic object, a significant part of the 
incident light can be scattered backward (reflected) what- 
ever the object be reflecting or transparent. It was re- 
cently demonstrated that an array of two-side-open sub- 
wavelength metallic slits effectively reflects light waves at 
the appropriate resonant conditions ^Li£ One may sup- 
pose that this is true also in the case of a single slit. 
In this section, we test whether a light wave can be 
resonantly reflected by a single two-side-open subwave- 
length metallic slit. To address this question, the energy 
flux in front of the slit is analyzed numerically for vari- 
ous regimes of the light scattering. In order to compare 
properties of the light reflection with that of the extraor- 
dinary (resonant) transmissioniii^i^, we consider the 
zeroth-order scattering of a time-harmonic wave of TM- 
polarized light by a slit in a perfectly conducting thick 
metal film placed in vacuum (Fig. The energy flux 
Si in front of the slit is compared with the fluxes Sn and 
Si 1 1 inside the slit and behind the slit, respectively. The 
amplitude and phase of the light wave at the slit entrance 
and exit are compared with that of a Fabry-Perot cav- 
ity. The electric E and magnetic H fields of the light are 
computed by using the traditional approach based on the 
Neerhoff and Mur solution of Maxwell's equationsiS^ 
For more details of the model, see the Appendix. 

According to the model, the electric E(x,z) and mag- 
netic H(x, z) fields in front of the slit (region I), inside the 
slit (region II) and behind the slit (region III) are deter- 
mined by the scalar fields U\{x, z), ^(x, z) and 1/3(2;, 2), 
respectively. The scalar fields are found by solving the 
Neerhoff and Mur integral equations. The magnetic field 
of the wave is assumed to be time harmonic and constant 
in the y direction: H{x, y, z, t) = U(x, z) exp(—iu>t)e y . In 
front of the slit, the field is decomposed into U\{x,z) — 
U l (x,z) + U r (x,z) + U d (x,z). The field W(x,z) rep- 
resents the incident field, which is assumed to be a 
plane wave of unit amplitude; U r (x,z) denotes the field 
that would be reflected if there were no slit in the film; 
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FIG. 1: Propagation of a continuous wave through a subwave- 
length nano-sized slit in a thick metal film. 



U d (x, z) describes the backward diffracted field due to 
the presence of the slit. The time averaged Poynting 
vector (energy flux) S of the electromagnetic field is cal- 
culated (in CGS units) as S = (c/16n)(Ex H* + E* x H). 
The reflection coefficient R — Sl d t is given by the nor- 
malized flux = S rd I ' S % integrated over the slit width 
2a at the slit entrance (z = b), where S rd is the z compo- 
nent of the backward scattered flux, and S % is the incident 
flux along the z direction. The flux S rd = S rd (U r ,U d ) 
is produced by the interference of the backward scat- 
tered fields U r (x,z) and U d (x,z). The transmission co- 
efficient T = Sf nt (b) is determined by the normalized 
flux S 3 = S 3 /S l integrated over the slit width at the slit 
exit (z — 0), where the flux z-component S 3 — S 3 (U 3 ) 
is produced by the forward scattered (transmitted) field 
U 3 (x,z). Notice that the definitions of the reflection R 
and transmission T coefficients are equivalent to the more 
convenient ones defined as the integrated reflected or 
transmitted flux divided by the integrated incident flux. 
In the following analysis, the reflection and transmission 
coefficients are compared to the fluxes Sf nt and Sf^t ob- 
tained by integrating the normalized fluxes S d — S d /S z 
and S% d = S ird /S\ respectively. 

We analyzed the backward scattering of light for a wide 
range of scattering conditions determined by values of the 
wavelength A, slit width 2a and film thickness b. As an 
example, the reflection coefficient R — S\ d t (b) as a func- 
tion of the film thickness b computed for the wavelength 
A = 800 nm and the slit width 2a = 25 nm is shown in 



Fig. 2(a) The transmission coefficient T = Sf nt (b) and 
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efficient correspond to reflection minima. In Fig. 2(a) 
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FIG. 2: (a) The reflection coefficient R — SX nt (6), the trans- 
mission coefficient T = Sf nt (b), and the integrated fluxes 
Sf nt (b) and S^f (6) as a function of the film thickness b com- 
puted for the wavelength A = 800 nm and the slit width 
2a = 25 nm. (b) The logarithm of the integrated flux 
Slnt ( a i b) as a function of the slit half- width a and film thick- 
ness 6. 



the integrated fluxes Sf nt (b) and 5-^(6) are presented in 
the figure for the comparison. We note the reflection res- 
onances of A/2 periodicity with the maxima R ma x ~ 2. 
In agreement with the previous resultsi&iiii 5 ^ one can 
see also the transmission resonances having the same pe- 
riod, and the peak heights T « 10 (T » A/27ra) at the 
resonances. It is worth to note the correlation between 
the positions of maxima and minima in the reflection and 
transmission. The resonance positions for the total reflec- 
tion are somewhat left-shifted with respect to the trans- 
mission resonances. The maxima of the transmission co- 



one can observe also many satellite peaks in reflection. 
For one broad minimum, it appears a local reflection 
maximum, which is characterized by a weak amplitude. 
The local maxima appear before 400, 800, and 1200 nm. 
The positions of the local maxima approximately cor- 
respond to the Sf nt maxima. To clarify a role of the 
fields U l , U r , U d and U 3 in the resonant backward scat- 
tering, we compared the integrated flux Sl^ d (U l , U r , U d ) 
with the fluxes Sj f (U d ) and Sf nt {U 3 ) = T. One can 



see from Fig. 2(a) that the flux produced in front of 



the slit by the interference of the incident field 17* (x, z) 
and the backward scattered fields U r (x,z) and U d (x,z) 
is practically undistinguishable from that generated by 
the backward diffracted field U d and forward scattered 
(transmitted) field U 3 . The integrated flux S\ r n d {a,b) 
as a function of the slit half-width a and film thick- 



ness b is shown in Fig. 2(b) We notice that the widths 
and shifts of the resonances increase with increasing the 
value a. Analysis of Fig. 2(a) indicates that the differ- 



ence between the integrated fluxes Sl d t (U r , U d ) = R and 
Sf nt (U 3 ) = T {T^Sf nt (U d )) appears due to the inter- 
ference of the backward diffracted field U (x, z) and the 
reflected field U r (x, z). 

The dispersion of the reflection coefficient R(X) — 

25 nm and the screen 
The inte- 



S\ d t {\) for the slit width 2a 



thickness b 
grated fluxes Sf nt 



351 nm is shown in Fig. 3(a) 



(A), S$(\) and S 3 nt (X)^ T(A) ver- 



sus the wavelength are shown in the figure for the com- 
parison. A very interesting behavior of the dispersion 
is that the coefficient R versus the wavelength A pos- 
sesses a Fabry-Perot-like dependence that is similar to the 
anomalous behavior of transmission T(A) reported in the 
studieo 10 i 13 i 14 i 15 i 19 . In agreement with the studiea 1 ' 2 ! 3 ' 15 , 
the height of the first (maximum) transmission peak is 
given by T « A 1 /2ira. The wavelengths corresponding 
to the resonant peaks A™ « 2b/m (m = 1, 2, 3, ... ) are 
in accordance with the resultsi 3 -. The high peak am- 
plitudes (enhancement), however, are different from the 
low magnitudes (attenuation) predicted in the studyi^, 
but compare well with the experimental and theoretical 
resu it s ij2 J 3ji4 i i5 difference is caused by the manner 

in which the Maxwell equations are solved. The studyi 3 - 
uses a simplified approach based on the matching the 
cavity modes expansion of the light wave inside the slit 
with the plain waves expansion above and below the slit 
using two boundary conditions, at z — and \z\ = b. 
Conversely, the Neerhoff and Mur method performs the 
matching with five boundary conditions, at z — > 0, z — > b, 
x — > a, x — > —a, and r — > oo. In contrast to the sharp 
Lorentzian-like transmission peaks, the slit forms very 
wide Fano-type reflection bands (see, Fig. 3(a) I. For one 
broad minimum in reflection, it appears also a local re- 
flection maximum, which is characterized by weak ampli- 
tude. At the near-to-resonance wavelengths of the trans- 
mission, the open aperture totally reflects the light. It is 
worth to note the correlation between the wavelengths for 
maxima and minima in the reflection R(X), transmission 
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FIG. 3: (a) The reflection coefficient R = S|jJ(A), the trans- 
mission coefficient T = Sf nt (A), and the integrated fluxes 
Sf n t(X) and Slnt (A) versus the wavelength A computed for the 
slit width 2a = 25 nm and the screen thickness 6 = 351 nm. 
(b) The real part of the normalized electric field x-component 
E X (U2) = E x (x,z) versus the normalized distance z/b in- 
side the slit cavity at x = 0, for the resonant wavelength 
A 1 = 800 nm, \l = 389 nm, A^ = 255 nm. 



T(A), and the flux Sf nt (Table [I]). The resonance wave- 
lengths for the main reflection maxima are red-shifted 
with respect to the transmission resonances. The wave- 
lengths of both the transmission and reflection (main 
and little) resonances are red-shifted with respect to 
the Fabry-Perot wavelengths A™ = 702 nm, 351 nm,. . . 
(A™ « 2b/ m, m = 1,2,3,...). To understand the 
physical mechanism of the resonant backward scattering, 
we also compared the integrated flux S\J^{U' l ,U r ,U d ) 
with the fluxes Sf nt {U d ) and Sf nt (U 3 ) = T. As can 
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TABLE I: The wavelengths for maxima and minima in the 
reflection R(\), transmission T(A), and the flux Sf nt . 



be seen from Fig. |3(a)| the integrated fluxes Sf nt (X), 
SlntW an d £f n t (A) = T(X) are practically undistin- 
guishable also in the A-domain (for the 6-domain, see 
Fig. 2(a) |. The difference between the integrated fluxes 
Sl* t \tr\U d ) = R and Sf nt (U 3 ) = T (T*Sf nt (U d )) is 
caused by the interference of the backward diffracted field 
U d (x, z) and the reflected field U r {x, z) in the energy flux 
S - (E x H* + E* x H). The wavelengths of the little 
maxima of the reflection R = R(U d , U r ) correspond ap- 
proximately to the high maxima S d nt . Therefore, the 
little maxima can be attributed to the interference of the 
reflected field U r with the dominant diffracted field U d . 
The red shifts and the asymmetrical shapes of the re- 
flection bands can be explained by a Fano analysis 8 of 
the scattering problem by distinguishing resonant and 
non-resonant interfering contributions to the reflection 
process. The resonant contribution is given by the field 
U d and the non-resonant one is attributed to the field 
U r . Other interesting interpretations of the shifts of res- 
onant wavelengths in the transmission spectra from the 
values 2b/ m can be found in the studie o 10 i 13 i 14 i 15 i 17 i 23 . It 
should be mentioned that the asymmetrical behavior of 
reflection was observed also in the case of a Fabry-Perot 
resonatoas^iS^. The conditions to achieve such an asym- 
metry are rested on the existence of dissipative loss in 
the resonator. There is no explicit loss in the present 
problem, but the dissipative loss can be substituted by 
radiative loss due to the diffraction by the slit. 

After the analysis of Fig. |2(a)| it is not surprising 
that the maxima of the transmission are accompanied 
by the minima of the reflection also in the A-domain (see, 
Fig. |3(a)] |. It should be noted in this connection that such 
a behavior of i?(A) and T(A) is similar to that observed in 
the case of excitation of the surface plasmons in an array 
of slit in a thin metal film. 18 In the study 18 , the minima 
in reflection spectra corresponding to the maxima in the 
transmission spectra were attributed to the redistribu- 
tion of the energy of diffracted evanescent order into the 
propagating order. In the case of a thick film, we explain 
such a behavior by another physical mechanism, the in- 
terference of the backward diffracted field U d (x,z) and 
the reflected field U r (x, z). It can be noted that the corre- 
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FIG. 4: (a) The phase distribution <p(x, z) of the electric field 
x-component E x (x,z) inside and outside the slit. The field 
component E x (x,z) is given by E x (U\U r ,U d ), E X (U 2 ) and 
E X {U 3 ) in the regions I, II and III, respectively, (b) The phase 
distribution <p(x, Z) at x = 0. The slit width 2a — 25 nm, the 
film thickness b = 351 nm and the wavelengths A = 800 nm. 



lation of positions of refl ect ion m inima and transmission 
maxima (see, Figs. 2(a) and 3(a) I are consistent with that 
predicted by the study 17 for TE-polarized light scattered 
by a grating of two-side-open slits in a thick metal film. 
However, the values of R(X) and T(A) are in contrast to 
the relation R(X) + T(A) = 1 given in the study 17 . The 
difference can be explained by the fact that we examined 
light scattering by an infinite screen using local defini- 
tions of R and T, while the study^ 7 . analyzed the global 
reflection and transmission by a grating of finite size. 

The dispersions Sj t (X), S**(\) and Sf nt (X) = T(X) 
shown in Fig. |3(a)| indicate the wave-cavity interaction 
behavior, which is similar to that in the case of a Fabry- 



FIG. 5: (a) The spatial distribution \S z (x,z)\ of the absolute 
value of the normalized energy flux along the z direction in- 
side and outside the slit. The distribution \S z (x,z)\ is shown 
in the logarithmical scale. The flux S z is given by the nor- 
malized fluxes S ird {U i ,U r 1 U d )/S i (U i ), S 2 (U 2 )/S i (U i ) and 
S 3 (U 3 )/S l (U i ) for the regions I, II and III, respectively, (b) 
The energy flux distribution S z {x, z) at x = 0. The slit width 
2a — 25 nm, the film thickness b = 351 nm and the wave- 
length A = 800 nm (solid line) and A = 882 nm (dotted line) 
corresponding to a transmission resonance and a little reflec- 
tion resonance, respectively. 



Perot resonator. The fluxes Sf nt (U d ), Sgf(U i ,U r ,U d ) 
and Sf nt (U 3 ) — T exhibit the Fabry-Perot-like maxima 
around the resonance wavelengths A™ » 2b /m. In or- 
der to understand the connection between the Fabry- 
Perot-like resonances and the total reflection, we com- 
puted the amplitude and phase distributions of the light 
wave at the resonant and near-resonant wavelengths in- 
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side and outside the slit cavity (see, Figs. 3(b) 0] and 
|SJ. At the resonance wavelengths, the intra-slit fields 
possess maximum am plitud es with Fabry-Perot-like spa- 
tial distributions (Fig. 3(b) \. However, in contrast to the 
Fabry-Perot-like modal distributions, the resonant con- 
figurations are characterized by antinodes of the electric 
field at each open aperture of the slit. Such a behavior 
is in agreement with the resultsi2*22i2i. It is interesting 
that at the slit entrance, the amplitudes E x of the res- 
onant field configuration possesses the Fabry-Perot-like 
phase shift on the value of n (Fig. @J. The integrated 
fluxes Sl r J(U\U r ,U d ), S d nt (U d ) and Sf nt (U 3 ), at the 
first resonant wavelength A*, exhibit enhancement by 
a fac tor A /27ra ps 10 with respect to the incident wave 
(Fig. 3(a) I. For the comparison, the normalized resonant 
fluxes S% d {U\ U r , U d ) and S 3 (U 3 ) in the near-field zone 
(z —2a) are about 5 times greater than the incident 
wave (see, Fig. It should be stressed that the res- 
onantly enhanced intra-cavity intensity S 2 (U 2 ) is about 
10 times higher than the resonant fluxes Sll d (U l 1 U r , U d ) 
and S 3 (U 3 ) localized in the near field zone in front of 
the slit and behind the slit, respectively (Fig. |3J|. The 
interference of the incident U l {x,z) wave and the back- 
ward scattered fields U r (x, z) and U d (x, z), at the reso- 
nant wavelengths A™ ~ 26/rn, produces in the near-field 
diffraction zone a strongly localized wave whose normal- 
ized flux S^ d (U\ U r , U d ) is A/27ra «10-10 3 times greater 
than the incident wave, but about one order of magnitude 
smaller than the resonant intra-cavity intensity. 

In our model we considered an incident wave with TM 
polarization. According to the theory of waveguides, the 
vectorial wave equations for this polarization can be re- 
duced to one scalar equation describing the magnetic field 
H of TM modes. The electric component E of these 
modes is found using the field H and Maxwell's equa- 
tions. The TM scalar equation for the component H 
is decoupled from the similar scalar equation describing 
the field E of TE (transverse electric) modes. Hence, 
the formalism works analogously for TE polarization ex- 
changing the E and H fields. 



terference of incident and resonantly backward-scattered 
light produces in the near-field diffraction zone a strongly 
localized wave whose intensity is greater than the in- 
cident wave by a factor A/27ra «10-10 3 and about one 
order of magnitude smaller than the intra-cavity inten- 
sity. The correlation between the amplitude and phase 
distributions of light waves inside and outside the slit 
was also investigated. The slit cavity was compared with 
a Fabry-Perot resonator. We showed that the ampli- 
tude and phase of the resonant wave at the slit entrance 
and exit are different from that of a Fabry-Perot cavity. 
The physical mechanism responsible for the total reflec- 
tion is the interference of the backward diffracted res- 
onant field U d (x, z) and the reflected non-resonant field 
U r (x, z) in the energy flux at the near-to-resonance wave- 
lengths (Fano-type effect). The wavelength-selective to- 
tal reflection of light by two-side-open metal slits may 
find application in many kinds of sensors and actuators. 
The (10-10 3 )-times and (10 2 -10 4 )-times enhancement of 
the light intensity in front of the slit and inside the slit 
can be used in reflective nanooptics and in intra-cavity 
spectroscopy of single atoms. We believe that the pre- 
sented results gain insight into the physics of resonant 
scattering of light by subwavelength nano-slits in metal 
films. 
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III. SUMMARY AND CONCLUSION 

In the present paper, the backward scattering of TM- 
polarizcd light by a two-side-open subwavelength slit in 
a metal film has been analyzed. We predict that the re- 
flection coefficient versus wavelength possesses a Fabry- 
Perot-like dependence that is similar to the anomalous 
behavior of transmission. The open slit totally reflects 
the light at the near-to-resonance wavelengths. The res- 
onance wavelengths for the total reflection are some- 
what red-shifted with respect to the transmission res- 
onances. The wavelengths of both the reflection and 
transmission resonances are red-shifted with respect to 
the Fabry-Perot wavelengths. The sharp resonant max- 
ima of transmission are accompanied by the wide minima 
of the reflection. In addition, we showed that the in- 



Appendix 

We briefly describe the Neerhoff and Mur modeli* 3 . 
of the scattering of a plane continuous wave by a sub- 
wavelength slit of width 2a in a perfectly conducting 
metal screen of thickness b. The slit is illuminated by 
a normally incident plane wave under TM polarization 
(magnetic- field vector parallel to the slit), as shown in 
Fig. 2] The magnetic field of the wave is assumed to be 
time harmonic and constant in the y direction: 

H(x,y,z,t) = U(x,z)e- lut e y . (1) 

The electric field of the wave is found from the scalar 
field U(x, z) using Maxwell's equations. The restrictions 
in Eq. Q reduce the diffraction problem to one involving 
a single scalar field U(x, z) in two dimensions. The field 



7 



is represented by Uj(x,z) (j— 1,2,3 in region I, II and 
III, respectively), and satisfies the Helmholtz equation: 
(V 2 + kfiUj = 0, where j = 1, 2, 3. In region I, the field 
Ui(x, z) is decomposed into three components: 

Ui(x,z) = U i (x,z) + U r {x,z) + U d (x,z), (2) 

each of which satisfies the Helmholtz equation. U % rep- 
resents the incident field: 

U*(x,z) =e~ lklZ . (3) 

U r denotes the reflected field without a slit: 

U r (x,z) = W(x,2b-z). (4) 

Finally, U d describes the diffracted field in region I due 
to the presence of the slit. With the above set of equa- 
tions and standard boundary conditions for a perfectly 
conducting screen, a unique solution exists for the scat- 
tering problem. The solution is found by using the Green 
function formalism. 

The magnetic H(x, z, t) fields in regions I, II, and III 
are given by 

H (x, z) = exp(— ik\z) + exp(— ik\{2b — z)) 

N 

ia t\ 

nT 2 



Y, Hh{kiJ(x - Xj f + (z- bf){DU b ) h (5) 



H 3 {x,z) 



N 



3=1 



xH, 



(i) 



- Xj) 2 + z 

a, J = 1,2, 



i (X) is the Hankel function; H l 



(7) 

,N; N > 
- H* • e,„ 



where Xj = 2a(j - 1/2) JN 
2a/ z- Hf 

i = 1,2,3; 7 m = [kl - (m7r/2a) 2 ] 1/2 . The coefficients 
(DUo)j are found by solving numerically four coupled 
integral equations. For more details on the model and 
the numerical solution of the Neerhoff and Mur coupled 
integral equations, see the references^. 



3=1 



H 2 (x,z) 



1 _ e V*tN 



JY 



2Nyn4 



2N v r k% 



N N 



3=1 

N 



3 = 1 



2N 



m=l ^ a 



3=1 



N 

X COS 
3=1 



mir(xj + a) 



2a 



(DU )j - - cos 



rmr{x + a) 



m— 1 



N 

E 

3=1 



COS 



mir(xj + a) 



rn— 1 



2a 

yi\\z—V\ 



N 



X COS ■ 



rmr(x + a) <r— » mn(xj + a) 



E 



cos ■ 



2a ^ 2a 

3=1 



-W)j 



1 00 _1_ 

- V cos(m7r^)e^l z - b l 



N 

N 



2a 

mir(xj + a) 



m— 1 

r 

cos ± (U b )j, 



3 = 1 



2a 



(6) 
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